Criticality in inhomogeneous magnetic systems: Application to quantum ferromagnets 
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We consider a (/)''-theory with a position-dependent distance from the critical point. One real- 
ization of this model is a classical ferromagnet subject to non-uniform mechanical stress. We find 
a sharp phase transition where the envelope of the local magnetization vanishes uniformly. The 
first-order transition in a quantum ferromagnet also remains sharp. The universal mechanism lead- 
ing to a tricritical point in an itinerant quantum ferromagnet is suppressed, and in principle one 
can recover a quantum critical point with mean-field exponents. Observable consequences of these 
results are discussed. 
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In standard phase transitions, such as the paramagnet- 
ferromagnet transition, or the liquid-gas transition, a ho- 
mogeneous order parameter (OP; the magnetization in a 
magnet, or the density difference in a fluid) goes to zero 
as one crosses from the ordered phase into the disordered 
one. The OP may vanish continuously, as in the case of a 
magnet where the transition is second order, or discontin- 
uously, as in the case of a fluid where the transition is first 
order except at the critical point. An external fleld may 
preclude a homogeneous OP. This happens for a fluid in a 
gravitational field, which produces a position-dependent 
density profile [l| and in some sense destroys the critical 
point (see below). Due to the weakness of gravity, this 
is a very small effect. This raises the question whether 
qualitatively similar, and maybe quantitatively larger, ef- 
fects can be achieved in other systems if an external field 
induces an inhomogeneous OP. 

We consider one such example, namely, a ferromag- 
net subject to mechanical stress. We will first discuss 
a classical Heisenberg magnet, and later generalize to 
quantum ferromagnets (FMs). Consider a metallic FM 
in the shape of a circular disk that is bent in the di- 
rection perpendicular to the disk plane. This leads to 
a position dependent mass density [2] and hence, in a 
metal, to a position dependent electron density and an 
inhomogeneous chemical potential /i. The FM transi- 
tion is described by a 0*-theory Q , and within a Stoner 
model the inhomogeneous /i leads to a spatially depen- 
dent distance from criticality. Naively, this means that 
the system can be tuned to criticality only at special po- 
sitions within the sample, not everywhere at the same 
time. One might thus expect the transition to become 
smeared. This is what appears to be found for the liquid- 
gas transition in a gravitational field which does, 
however, present a physically different situation For 
the FM case we find that the transition remains sharp in a 
well-defined sense with mean-field critical behavior, even 
though the magnetization Ai'(x) is position dependent 
and hence "smeared" in some sense [Sl, '9'| . One might ex- 



pect that M{x) is essentially restricted to a surface layer 
of fixed width, so that the dimensionality of the system 
is effectively reduced by one. This is not the case; we 
find (for a particular model of a sample of linear dimen- 
sion L) that Af(x) is essentially nonzero in a region of 
width L^/^, so the support of the magnetization diverges 
in the thermodynamic limit, L oo, albeit more slowly 
than L. This leads to unusual critical exponents for some 
spatially averaged observables. 

Consider a flat circular disk sample of a metallic FM 
with thickness L (in the x-direction) and a radius that 
is some flxed multiple of L. A distortion of the disk in 
x-direction from a flat shape into a paraboloid leads to 
a strain tensor whose trace is a linear function of a; 0], 
and hence to an electron density n{x) = no + const, x x. 
(This distortion must be achieved by bending, not, e.g., 
by grinding.) Within a Stoner model, the distance r from 
criticality depends linearly on the density of states, and 
hence on the cube root of n{x). We consider a more 
general model where r varies as a power of x: 



r{x) = ro + 2{xlLf 



(1) 



Note that we take the prefactor of the x-dependent term 
to be of 0(1) in order to demonstrate the qualitative ef- 
fects of such a term. For a real bent plate the prefactor 
will be smaller, and we will give a semi-quantitative dis- 
cussion below. Also note that r(L) is bounded as L — > oo. 
This reflects a bending displacement proportional to L 
and ensures that a meaningful thermodynamic limit can 
be taken. Our model action is a (^"^-theory with a spa- 
tially dependent mass given by r{x), 



S = 



dx 



^(V</.(=r))^ 



{cf,\x)y 



(2) 

The integration extends over a volume V oc L'^, and c 
and u are constants. We emphasize that this model is 
rather general, and a magnet under stress is only one 
possible realization. We flrst treat this problem in a 
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saddle-point approximation and look for solutions of the 
form 4>{x) = (0, 0, M{x)), with M{x) the inhomogeneous 
magnetization. The saddle-point equation then reads 

cM"{x) = r{x) M{x) +uM^{x) - H, (3) 

where we have added a magnetic field H. The physical 
solution of the ODE ^ must obey the boundary condi- 
tions M'(0) = M'{L) = 0. 

While this ODE would be difficult to solve in closed 
form, we can obtain a great deal of information from 
asymptotic solutions and scaling considerations. For x <C 
L one can neglect the x-dependence of r{x) and finds a 
solution in terms of the Jacobi integral sn(a;); for x ~ L 
one can drop the term and finds a solution in terms 
of Airy functions. A smeared transition would imply a 
nonzero magnetization for all parameter values This 
is physically not possible: for rg > 0, r{x) is positive 
definite and the physical solution must be M{x) = 0. 
We thus expect a sharp phase transition in the following 
sense: there exists a value rg of vq such that the envelope 
of the magnetization vanishes uniformly as tq — > rg at 

= 0, or as -ff — > at ro = i'q. 

There are two explicit length scales in this problem: 
the zero of r(x), xo — L(— ro/2)^/", and the bare cor- 
relation length ^0 = c/ro. One expects the phase 
transition to occur when xq = £,o (apart from a factor of 
0(1)). This condition leads to rg = _2/£2"/("+2)^ ^itji 
£ = i/Vc a dimensionless system size. For ro = rg the 
zero of r{x) is xg = L/^^/f""'"^'. Now define y = x/xq, 
^{y) = M{x), and h = Hl^"/^"+^\ ^ obeys 

(y) = Piv, ro^'"/("+') ) M(y) + (y) - K (4a) 

where 

p{y,z)^z + 2y^. (4b) 
The solution of Eq. (|^a|) determines M{x) via the relation 
M(a;) =r"/("+2)^(a;/a;g;ro^2n/(n+2)^^^3n/(n+2)^ ^ 

(5) 

where we show the dependence of jj. on ro and H . 

Now consider the thermodynamic limit, € — > oo. Since 
r{x) is bounded for all x (see Eq. ([1])), we expect physical 
quantities in this limit to be independent of I. From Eq. 
([5]) , this requirement yields for the local magnetization a 
power-law prefactor times an envelope function, 

M(x;ro,i/ = 0) = r^' 5;;j(a;/xg), (6a) 
M(x;ro-0,H) = i/^/^ ^^.(x/xg). (6b) 

Similarly, for the envelope susceptibility x{^) = 
{dM{x) / dH)H=o we have 

X(x; ro) = r^'^g^ix/x^), (6c) 

where g^, g^j, and g^^ are scaling functions. We can 
further define a local specific heat C{x) = d'^f{x)/drQ, 
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FIG. 1: Tlie a;-dependent magnetization (solid line) for ro = 
-1, n = 1, L = 1, c = 1/25, u = 1, H ^ 0. The dashed and 
dotted lines represent asymptotic analytic solutions for x <^ L 
and X ft: L, respectively. The inset shows the magnetization 
together with the x-dependent mass function r{x). 

where f{x) is the free energy density which scales as 
AI'^(x). For the critical exponents /3, 6, 7, and a de- 
fined by M{x) cx rp, M{x) cx H^^^, x{^) oc i"q^ , and 
C{x) cx r^" this implies 

^=1/2 , ,5 = 3 , 7=1 , a = 0. (7a) 

Finally, we determine the exponents and 77. The mag- 
netization depends on rg only through the combination 
ro ("+2) ^ If identify the diverging length scale ^ 
that characterizes the phase transition with c^/^£i'^/i^+^) 
we have ro ~ or ^ ~ r^'^ with u = 1/2. Further- 
more, the inverse susceptibility determines the exponent 
77 via cx ro ^ = ^^^+1) -^^i^jj = q \Ye thus have 

, 77 = 0. (7b) 

1/ and 77 defined in this way are finite-size scaling expo- 
nents, and v does not represent the divergence of a co- 
herence length defined via the spatial decay of the two- 
point correlation function. The latter remains finite even 
at the transition, as it does in the case of the liquid-gas 
transition in a gravitational field Q . There thus is also a 
correlation length exponent v that is equal to zero. There 
is a further ambiguity within the framework of finite-size 
scaling: defining ^, for instance, as ^ = c^^^i would lead 
to different values of v and 77. The mean-field values 
given above result from what in some sense is the most 
natural choice for ^. However, all choices for ^ preserve 
the exponent relation i/(2 — 77) — 7. The Essam-Fisher 
relation, a + 2/3 + 'y = 2, is also fulfilled. 

We have corroborated and augmented these results by 
solving Eq. ([3]) numerically with the boundary condition 
M'(0) = 0, and M(0) chosen such that the mean-field 
free energy is minimized. Fig. [T] shows M(x) for the case 
of a linearly x-dependent mass, n — 1. M{x) is large in 
the region where r{x) < and small in the region where 
r{x) > 0, as one would expect. The asymptotic solutions 
mentioned above are also shown in Fig. [T] 

With increasing ro , the magnetization decreases for all 
X, see Fig. [2J and it vanishes uniformly when ro reaches 
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the critical value Tq. For tq > the physical solution 
of Eq. ^ is M{x) = 0. As ro — > rfj from below, the 
envelope of the magnetization vanishes as I^q — 
For the value M{x = 0) this is demonstrated in the inset 
in Fig. O In a magnetic field the magnetization vanishes 
uniformly as H^^^ for H at tq = fg, see Fig. [3] 
These results show that there is a sharp phase transition: 
the envelope of the magnetization vanishes uniformly as 
ro ~f Tq a.t H = 0, or as if ^ at ro = rg. The order- 
parameter critical exponents have mean-field values: /3 = 
1/2 and (5 = 3. The numerics suggest that the mean-field 
critical behavior also holds for finite L. We have found 
Eqs. ([6al I6bp to hold for values of L^/c as small as 16. 

We now consider spatially averaged observables (de- 
noted by an overbar) rather than local quantities. Con- 
sider an averaged magnetization M — (l/L) dx M{x). 
Since M{x) is essentially nonzero only on the interval 
X e [0, ig], the upper limit of the integral is essentially 
Xq oc Corrections to this approximation arc 

exponentially small. M obeys, instead of Eq. (O, 

M = r^/Af (ro^'"/("+'),i/^3"/("+')) ,(8a) 



fM{u,v) 



dy^i{y;u,v) 



(8b) 



Demanding again that observables are independent of £ 
for £ oo, this leads to exponents p = {n + 2)/2n and 
6 = 3n/{n + 2). Analogous considerations for the other 
observables we have considered yield the following set of 
exponents for spatially averaged quantities, 



2n 



3n 



1 



1 = 



— . (9a) 
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From the averaged susceptibility x = {DM /dH)H=o we 
find the exponent fj, and v is unchanged since M depends 
on the same combination of tq and £ as M{x), 



(9b) 
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FIG. 2: The a;-dependent magnetization for ro = —1 (1), 
ro = —0.8 (2), and ro = —0.6 (3). All other parameter values 
are the same as in Fig. [1] which lead to rg = —0.55213. The 
inset shows a log-log plot of M{x = 0) vs. ro — ro, with the 
solid line representing a power law with an exponent of 1/2. 
See the text for additional information. 
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FIG. 3: The x-dependent magnetization for ro = rg and H — 
0.003 (1), H = 0.0005 (2), and H = 0.00005 (3). All other 
parameter values are the same as in Fig. [1] The inset shows a 
log-log plot of M{x = 0) vs H , with the solid line representing 
a power law with an exponent of 1/3. 



These exponents for the averaged quantities satisfy again 
the relations a + 2(3 + j = 2 and iy{2 — fy) = 7. 

The exponent values depend on how the averaged 
quantities are defined, and hence on what exactly is 
being measured. For instance, if one defined M as 
(l/c^/^L^/^) dx M{x) to account for the fact that the 
magnetization is essentially nonzero only for x < Xq (x 
£^^^, one would find mean-field values for all exponents. 

We now apply these results to quantum FMs. It is ob- 
served that the transition in itinerant FMs at sufficiently 
low temperatures T is always first order |1(|. This has 
been explained in terms of fluctuation effects due to the 
coupling of the OP fluctuations to soft particle-hole ex- 
citations in metallic FMs. As a result, the free energy in 
a mean-field approximation has the form [11 1 



F = - Af 
2 
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©(Af*^), (10) 



with t; > 0. The v-term is negative, which leads to a 
first-order transition at T = 0, and to a tricritical point 
at a temperature Ttc = exp(— w/2t)) [l^ . 

A position dependent chemical potential in the regime 
where the transition is first order has two effects: (1) 
It lowers the tricritical temperature Ttc, and a suffi- 
ciently strong space dependence restores a quantum criti- 
cal point; (2) it leads to an inhomogeneous magnetization 
as in the case of classical magnets discussed above. 

To quantify the first effect we estimate the prefactor 
of the space dependent term in /i. For a displacement 
at the edge of the plate equal to 0.01 L (conservatively; 
more severe bending may be possible) one finds a den- 
sity variation of about 0.01 rtg, and, instead of Eq. ([1]), 
r[x) — ro + 2px/ Lwii\\ p — 0(0.01). Suppressing factors 
of 0(1), the chemical potential is p,{x) — e-p{l + px/L). 
The effects are smaller by a factor of about 10^ compared 
to those shown in our figures, but still large compared to 
those of gravity on a fluid. To estimate the effect on 
Ttc we recall that the ii-term in the free energy is due 
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FIG. 4: M{x = 0) vs. ro (solid line, left scale) as obtained 
from the solution of Eq. [TT]for parameter values L — 1, c — 
1/25, u = -1, w = 1, which lead to = -0.39216. Also 
shown is the free energy density / (dashed line, right scale) . 



to a soft propagator P = l/{oj + kvp) [13], with oj (k) 
the frequency (wave number). A position dependent /j, 
causes P to acquire a mass or cutoff frequency iVc pro- 
portional to the prefactor of the a;-dependency in fj,{x), 
ujc ~ p ep / hLkp . Ttc is lowered by roughly huc/kB- With 
p as above, and typical values for all other quantities, one 
finds a very small suppression i5Ttc = O(10~^K). It is 
interesting, however, that a sufficiently strong spatial de- 
pendence of p will eventually destroy the mechanism that 
leads to a first-order quantum phase transition, as does 
quenched disorder [l^. 

To study the second effect we thus assume that the 
transition is still first order at low temperature. In this 
case the logarithm in Eq. pO|) can be expanded, and the 
Landau free energy is adequately represented by a power 
series with a negative M'^ term and a positive term. 
Generalizing to a non-homogeneous situation we then ob- 
tain the following ODE for the OP in the quantum case, 

cM"{x) = r{x)M{x)+uM^{x)+wM'''{x), (11) 

which replaces Eq. ([3]). Here w > 0, w < 0, and r(x) 
IS given by Eq. H]). The resulting M{x) curves are very 
similar to those in Fig. [21 but now M{x) goes to zero 
discontinuously at a critical value Tq of tq, see Fig. |31 

We now consider the validity of our mean-field treat- 
ment. If the quantum critical point is restored by a (hy- 
pothetical) large position dependence of p, the quantum 
critical behavior will be mean-field like since the spatial 
dependence of p suppresses the mechanism that causes 
a first-order transition in clean quantum FMs, and non- 
mean- field critical behavior in disordered ones [l3| . For 
the classical transition it is possible that the mean-field 
critical behavior will also be exact. This is because the 
two-point correlation function is always finite, and hence 
there are no divergencies in simple perturbation theory. 
However, the divergence of the envelope susceptibility, 
Eq. (fSc]) . shows that there are fluctuations that may in- 
fluence the critical behavior. Fluctuations of the elas- 
tic deformation should also be considered 14 , 1^ ■ These 
points require additional investigation. 



In summary, we have considered a model for a metallic 
FM with a position dependent electron density or chemi- 
cal potential. This can be realized by mechanically stress- 
ing the sample. The phase transition remains sharp even 
though the electron density is not homogeneous, and we 
have given critical exponents for both local and spatially 
averaged observables. In the quantum case, the tricriti- 
cal temperature is lowered, although for realistic stresses 
this is a small effect. If a stronger position dependence 
can be realized (e.g., by means of optical lattices), this 
will eventually suppress the tricritical point, restoring a 
quantum critical point with mean-field critical behavior. 
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